IDENTIFICATION AND USE OF INFORMATIVE GENETIC SEQUENCES 



CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] The present application claims the priority benefit of p e nding U.S. Provisional Patent Application 
No. 60/441,745 fil e d January 23 , 2003 and pending U.S. Provisional Patent Application No. 60/441,806 
filed January 23, 2003. The complete disclosure of both-that applications is hereby incorporated herein 
by reference in its entirety . The pr e s e nt application also h e r e by incorporate s herein by reference th e 

entire disclosure of p e nding utility pat e nt application 10/ , titled M e thod and Syotcm for 

Id e nti f ying Biological Entiti es in Biological and Environmental Sampl e s, fil e d January 23, 200 1. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT 

[0002] Embodiments of the invention may have been conceived or first actually reduced to practice in 

the performance of work under the following Government contracts: . ,- As a 

result, the Government may have certain rights in those inventions. 

REFERENCE TO DOCUMENTS CO-FILED ON CD-ROM 

[00031 A total of two identical CD-ROM disks ( labeled "Copy 1" and "Copy Tl are submitted he rewith 
each containing the following electronic files. The CD-ROMs were created on October 12. 2006 and the 
sizes of each file are listed parenthetically as follows. Each CD-ROM contains one file of sequen ce data 
entitled SEQUENCE LlSTING.txt (4,015 KB). All electronic files on these CD-ROM disks are herein 
incorporated bv reference in their entirety. 

BACKGROUND 

fQQQ3}f00041 Field of the Invention . Embodiments of the invention relate to the identification of 
genomic sequences that are informative of the biological characteristics (e.g., presence, abundance, 
virulence, genetic modification) of a sample, along with systems and methods of using such sequences for 
gathering information on one or more sets of organisms present in the sample. Specific embodiments 
relate to microbial organisms. 
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fQfrWrOOl 11 Embodiments of the invention include methods for identifying oligonucleotide sequences 
unique to a set of organisms. In those methods, after identifying unique genomic sequences (typically, 
but not necessarily, as described above), unique genomic sequences are divided into target-length 
oligonucleotide sequences; preferable a number of sequences starting with the beginning of the genomic 
sequence and progressing one nucleotide at a time. The oligonucleotide sequences are properly formatted 
and searched to a selected database using a similarity search engine such as BLAST. The results are 
parsed for those oligonucleotides showing uniqueness to the set. A computer program product can 
implement this method in modules similar to those above. 

[00 111 TOO 121 The invention includes methods for inferring both unique genomic sequences and unique 
oligonucleotide sequences. In those methods, genomic data characteristic of a first set of organisms is 
obtained. The data is formatted into at least one query-length sequence, each query-length sequence 
being of a format compatible with a similarity search engine. A selected genomic database, e.g., an NCBI 
database, is searhed using the query and the similarity search engine such as BLAST. The results are 
parsed for those sequences not associated with the first set of organisms, but showing similarity beyond a 
threshold. 

BRIEF DESCRIPTION OF THE DRAWINGS 

fQ£4^1[00131 Figure 1 is a flowchart describing, in conjunction with portions of the written description, 
methods of the present invention. 

f094»tr00141 F igure 2 is a notional illustration of a target oligonucleotide window in the context of a 
unique genomic sequence of the present invention. 

fOO444r00151 Figure 3 is a notional illustration of differential hybridization of several unique 
oligonucleotides of the present invention between two strains of E. coli. 

fQQ4£}rQ016 ] Figure 4 is an illustration of a decision tree of the present invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

£004-641 00 171 As required, detailed embodiments of the present invention are disclosed herein. 
However, it is to be understood that the disclosed embodiments are merely exemplary of the invention 
that may be embodied in various and alternative forms. The figures are not necessarily to scale, and some 
features may be exaggerated or minimized to show details of particular components. Therefore, specific 
structural and functional details disclosed herein are not to be interpreted as limiting, but merely as a basis 
for the claims and as a representative basis for teaching one skilled in the art to variously employ the 
present invention. 

(004^00181 For purposes of the invention disclosure, the term "primer" includes a short pre-existing 
polynucleotide chain to which new nucleotides can be added by DNA or RNA polymerase. 

fOO4Sir00191 The term "randomly amplifying" includes increasing the copy number of a segment of 
nucleic acid in vitro using random primers. 

f004CT00201 "Amplicon" refers to DNA that has been manufactured utilizing a polymerase chain 
reaction (PCR) where a set of single stranded primers is used to direct the amplification of a single 
species of DNA. 

fQQgffl [002 1 1 "Biological entity" describes a biological element, cellular component , or organism that 
exists as a particular and discrete unit. This includes, but is not limited to gene, transgene, oncogene, 
allele, protein, DNA, RNA, mitochondria, pathogenic trait, vector, plasmid, clone, Acytota, prokaryotes, 
eukaryotes, Protista, Fungi, Plantae, Animalia and Monera, or any mixture thereof. For simplicity, 
"organism" is used interchangeably herein with "biological entity," and includes sub-organism entities. 
"Set of organisms" includes sets of one biological entity. 

f0021i r00221 A "sample" may be from any source, and can be a gas, a fluid, a solid, a biological 
sample, an environmental sample, or any mixture thereof. 

fQQ3?}f0Q231 "Nucleic acids" means RNA and/or DNA, and may include unnatural bases. 
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[00231 [00241 The "unique oligonucleotide sequence" generally identifies a nucleic acid sequence for 
which the sequence is known and determined unique to a set of organisms. In some embodiments, unique 
oligonucleotide sequences are more than 30 nucleotides in length. 

fflQ24}|-()Q25] The terms "unique genomic sequence" and "unique sequence" are interchangeable in the 
invention and refer generally to a sequence of nucleic acids that are specific to a set of organisms. 

r0025ir00261 In the literature there exist at least two confusing nomenclature systems for referring to 
hybridization partners. Both use common terms: "probes" and "targets." For the purpose of this 
disclosure, a "target" is the known oligonucleotide sequence (preferably set-unique), whereas a "probe" is 
the nucleic acid sample whose characteristic(s) (e.g., identity, abundance, virulence) is being detected. 
"Probe" includes any single stranded nucleic acid sequence, molecule, genomic sequence, or amplicon 
which typically is labeled. Probes can hybridize to a target if sufficient complementarities exist. Note 
that labeling can be implemented at various stages in either the probe or target or both, as known to those 
skilled in the art. 

fOQ2frH00271 The terms "microarray" and "array" are interchangeable and include a set of miniaturized 
chemical or biological reaction areas that may also be used to test DNA, DNA fragments, RNA, 
antibodies, or proteins. 

fOOCTr00281 A "labeled" or "detectable" nucleic acid is a nucleic acid that can be detected. The term 
"detection" refers to a method where analysis or viewing of the detectable nucleic acid is possible visually 
or with the aid of a device, including, but not limited to microscopes, fluorescent activated cell sorter 
(FACS) devices, spectrophotometers, scintillation counters, and fluorometers, devices using mass 
spectrometry, devices using radio isotopes. 

fQQ2&p0291 "Hybridized" means having formed a sufficient number of base pairs to form a nucleic 
acid that is at least partly double-stranded under the conditions of detection. The term "hybridization" 
refers to the process by which two complementary strands of nucleic acids combine to form double- 
stranded molecules. 
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f 00291 [00301 The term "complementarity" refers to a property conferred by the base sequence of a 
single strand of DNA or RNA that may form a hybrid or double stranded DNA:DNA, RNA:RNA or 
DNA:RNA through hydrogen bonding between base pairs on the respective strands. Adenine (A) usually 
complements thymine (T) or uracil (U), while guanine (G) usually complements cytosine (C). 

fOQ3fflfOQ311 The identification and characterization of microorganisms present in the environment has 
historically been accomplished by exploiting a variety of biological, immunological, biochemical, and 
genetic differences between organisms. Nucleic acid-based diagnostic methods have been developed that 
are specific for a single organism or small sets of organisms. PCR-based assays are typically performed 
by designing oligonucleotide primers that amplify organism-specific fragments of DNA. These 
fragments are subsequently detected by methods such as gel-electrophoresis, real-time PCR, or 
hybridization to either a membrane or microarray. A limitation of these existing assays is that although a 
positive result is informative for a specific organism or organism set, a negative result typically provides 
little or no information about the organism(s) under investigation. 

fQQ3-Hf 00321 Though it is possible to multiplex primers for numerous amplifications in a PCR for the 
concurrent identification of a variety of organisms, it is non-trivial to design compatible multiple primer 
pair sets that function in a single amplification reaction. Thus, the number of microorganisms that can be 
detected or otherwise characterized concurrently with this type of multiplex reaction is relatively small. 
Techniques such as real-time PCR and quantitative PCR are limited in the number of primer sets that can 
be used in a single amplification reaction and in the number of fluorescent molecules available for 
labeling DNA molecules and detection. 

fQO3£|f00331 In a method for detecting and distinguishing between various species and strains of 
viruses, viral RNA is reverse transcribed from semi-random primers, amplified by specific primers and 
then labeled with fluorescent nucleotides in a non-amplifying reaction. The labeled nucleic acids are then 
hybridized to microarrays that have been spotted with virus and strain-specific oligonucleotides that 
represent the entire genomes of these organisms. The resulting hybridization pattern discriminates 
between viruses represented on the array. However this approach is not directly translatable to fungi and 
bacteria. The relatively large size (3-5 million bases) and complexity of bacterial and fungal genomes, as 
compared to most viral genomes, represents an obstacle in the ability to identify oligonucleotides that are 
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species and strain specific. In addition, it is not feasible to synthesize and spot every possible 
oligonucleotide for every microbial genome onto a microarray. 

fQO^r00341 Bioinformatic tools such as BLAST, are intended to identify similarities between 
sequences. While similarities between the sequences of organisms are useful in some types of analysis, 
the differences between genomes can be useful in the identification and characterization of 
microorganisms. Unfortunately, bacterial and fungal genomes are so large that it is resource-intensive to 
subtract common sequences in order to identify unique sequences from all known genomes. Frequently 
only small fragments of genomic sequences have been identified as unique are available for identification 
of an organism. The increasing number of genomes that have been, or will soon be, sequenced is one 
incentive for identifying large fragments of known genomes. 

fQQMir00351 Current nucleotide-based methods of identifying organisms typically rely on primer- 
requiring multiples PCR methods, or oligonucleotide microarrays that utilize the limited amount of 
ribosomal genes (approximately 1% of the genome), or costly shotgun approaches directed at entire 
genomes. Ribosomal DNA is highly conserved across higher level sets - therefore the amount of unique 
sequence is limited. 

fQQ££H"00361 Approximately 300 microbial genomes have been partially or completely sequenced to 
date. In spite of this wealth of information, existing methods of detection and characterization of 
microbes are limited by the availability of unique sequence information from the genomes of these 
organisms. Typically, only small fragments of genomic sequences are identified as unique. 

fQ93&irQ0371 Current DNA amplification approaches to identify microorganisms are limited in terms 
of the number of sequences that can be identified concurrently. In vitro, two methods are used to 
multiplex, or identify multiple sequences concurrently. Both are limited by the challenge of generating 
specific primer pair sets that work well together in a single reaction mixture. One method for assessing 
which amplicons are produced in a multiplex PCR is to run the amplification product on a gel and to 
discriminate the various amplicons based on molecular size. The number of size indicative bands that can 
be resolved on the gel is a limiting factor for this approach. 
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fQ(&74f003S1 Another method, real time PCR, uses different fluorescent tags to identify specific 
amplicons in a multiplex PCR. The number of amplicons that can be resolved using this approach is 
limited by the number of different fluorescent tags available for probes used in the reaction. Current 
limitations for fluorescent resonance energy transfer methods, such a Taqman TAQM.AN© PCR 
technology and molecular beacons are about four amplicons for a single reaction. Thus the limitations are 
at least two-fold: the first is a compatibility issue regarding the use of multiple sets of unique primers; and 
the second is resolution of the amplified products. 

Unique genomic sequences in a set of organisms (which set can consist of one or more 
biological entities) may include both coding and non-coding sequences. Coding sequences are sequences 
that are further processed into proteins or polypeptides, typically performing a single function. These 
sequences are frequently conserved across genus and species. Conserved coding sequences can include 
genes that code for enzymatic elements, structural elements, virulence factors or developmental specific 
functions and processes. 

ffiO3£}[004()) Non-coding sequences are sequences that are not further processed and do not appear to 
possess a known function at this time. These sequences may be contained in a portion of the genome that 
contains unique coding sequences as well as between conserved coding sequences. Since non-coding 
sequences do not provide a known function, they are frequently overlooked as unimportant genomic 
material. These unique non-coding sequences can be used to identify an organism, just as unique coding 
sequences are used. Informative sequences can reflect a variety of features e.g. structural, functional, 
metabolic, virulence. Set-unique sequences can be coding or non-coding sequences. Set-unique 
sequences (coding or non-coding) can be inferred (see below) or found by searching through fully 
sequenced genomes. Partially sequenced genomes typically focus on coding sequences. Combinations of 
sequences that are not necessarily individually set-unique can also be informative. Sequences unique to 
sets above the species or strain level can be bio-informative, e.g., used in analyzing a sample for 
information about the organisms in the sample as described below. 

fOQ40j[004 1 1 Embodiments of the present invention include methods and systems for the identification 
of genomic sequences that are informative of the biological characteristics (e.g., presence, abundance, 
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virulence, genetic modification) of a sample. Referring to Figure 1, a illustrative method 100 of the 
present invention is shown. 

[004 11 [00421 In the illustrated embodiment, a subset of the genomic data 105 of the organism under 
investigation is obtained. The subset 105 can be obtained from known genomic data source 10, e.g., 
UniGene, GenBank GENBANK® genomic data source , European Molecular Biological Laboratory 
(EMBL), among other sources. Genomic data can also be obtained as sequence derived from in vivo or in 
vitro experiments 20 such as PCR and enzymatic digestion. A preferred subset of genomic data is the 
entire genomic sequence of an organism. 

[0042] [00431 In some embodiments, the obtained genomic data is preprocessed 110. Each aspect of 
preprocessing can be performed as needed or desired. 

[00431 [00441 In preferred embodiments, if necessary, the genomic data subset is converted from its 
native format, e.g., standard GeaBm kGENBANK annotated format, to a format compatible with 
subsequent steps. In some embodiments, where GettBmk GENBANK genomic data source annotated 
form at is used, the genomic data is converted to FASTA format to support a subsequent BLAST search. 

fOO44ir00451 P reprocessing can involve removing or masking portions of the genomic data that are 
judged not likely to have informative value. In preferred embodiments, these portions are removed, partly 
to make the subsequent search more efficient. This can include sequences known to be conserved with 
respect to the organism set under investigation (though some conserved sequences can be bioinformative, 
and sequences conserved at one tier in a taxonomy may be unique in another tier), repeats, inverted 
repeats, long terminal repeats, sequences otherwise known to be not favorable for hybridization. 

[Q0 4 51 TO0461 In preferred embodiments, genomic data is divided into query-length sequences 115. 
Some embodiments start with sequences of 1000 bases in length. The speed of the search is one factor in 
selecting the size of the initial query sequence. Subsequent iterations, described below, divide the initial 
query length further. 

fQQ46P0471 In further preprocessing, query-length genomic sequences were realigned with the 
genome from which they were generated in order to determine the start and stop point of each query 
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length sequence within the genome. Any annotations within the genome in the region containing the 
query length genomic sequence were transferred to the query length genomic sequence. Annotated 
regions include sequences known to have a specific biological function such as protein coding regions, 
biologically active RNA encoding regions, promoter and regulatory elements, spacing elements within 
operons, protein binding sites, and the like. 

£0047^0048! Note that if an entire genome was obtained, and no preprocessing performed, the query- 
length sequence 115 is the entire genome of the organism under investigation. In some embodiments, all 
the genomic data available for the organism under investigation is obtained, all preprocessing steps are 
completed, resulting in annotated query-length sequences of 1000 bases that do not include conserved 
sequences, repeats of various types, or sequences having characteristics that otherwise make them not 
amenable to subsequent steps. 

feQ48¥0Q491 In preferred embodiments, the query length sequence (preprocessed or not) is used as a 
query to a similarity search program 120, e.g., BLAST. The query is directed to a selected database 125 
of genome data. In some embodiments, the selected database is limited to organisms of the same general 
type as that under question, in order to increase search efficiency over what it would be were the search 
directed to a full database containing a broader variety of organisms. In some embodiments, the query is 
directed to the NCBI nr database. For example, if only microbial organisms were under investigation, the 
selected database 125 would be a database of microbial genomic data - broader databases including, for 
example, mammalian genomic data, would be avoided at this stage. In these circumstances, a subsequent 
search against the broader database is preferred in order to confirm the set-uniqueness of these initial 
results. In some embodiments, query-length sequence is removed from the selected database, while in 
other embodiments, results showing homology to the query itself are either ignored, or taken as 
confirmation of the validity of the query with respect to the organism under investigation. Table 1 lists an 
exemplary query for a strain of Escherichia coli. 



| SEQ ID NO:l >NC 000913 29 part354 Escherichia coli K12, complete genome. 
TATGTTTTAAATAACTAATTGGTCGGGTTAGTGCATCCGGCTTTCTTTATATTCGCCAGA 
AGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGG 
T C AGGGTTGGGCGGCAGATAT C AC AGAGGTAGAAAC CAC C ACAGGTGAAAAGAAAAATAC 
CAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACC 
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CTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGC 

TGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCG 

CAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCC 

CCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCC 

CAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTAACTATT 

CGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCA 

ATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCT 

TGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTT 

ATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATT 

CCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGT 

CAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGC 

ACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGT 

CAATAACGAAGGCGAGAGCACCATCACCAAC 



TABLE 1 



fQQ494r00501 Preferred embodiments parse 130 the similarity search program output 125 to identify 
sequences lacking significant similarity with other organisms in the selected database, e.g., unique 
genomic sequences 132. This is counter to the typical use of such search programs. For the purpose of 
this disclosure, "unique" or "uniqueness" is a function of thresholds, preferably controlled by the user, 
regarding identity, homology, score, expected (E) value and the length of the unique sequence under 
consideration. Identity, score, expected value are data returned in a typical BLAST search. In some 
embodiments, lacking significant similarity, e.g., "unique," means no BLAST hits or hits with a E-value 
less than le-5. 

f£Q5#}|"005n In practice, computational resources are finite, so the selected database may range from a 
database of all fully or partially known genomes to a narrower database such as known microbial 
genomes. Directing the initial query to a database of less than all available genome data, while 
computationally economical, can make it advisable to BLAST the candidate sequences (e.g., in preferred 
embodiments, those genetic sequence segments found to be unique) against the broader databases, e.g., 
the NCBI nr database, to detect homology with known genomes. 

£Oe#Hf0052J_Table 2 is an extract from an exemplary BLAST output for the sequence in Table 1. Note 
that significant alignments with E values less that le-47 are all E. coli. This confirms that the query 
sequence is sufficiently "unique" with respect to E. coli to be informative. 
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BLASTN 2.2.7 [Jan-02-2004] 

RID: 1074618061 -2346 8- 21421 92 4 927. BLASTQ4 
Query= (1000 letters) 



Database: All GenBank+EMBL+DDBJ+PDB sequences (but no EST, STS, 

GSS, or phase 0, 1 or 2 HTGS sequences) 2,017,250 sequences; 9,771,119,756 total letters 
Distribution of 64 Blast Hits on the Query Sequence 



Sequences producing significant alignments: bits) 

SBQ ID NO: 2 gi | 1787665 | gb | AE000237 . 1 [AE000237 Escherichia coli K12 MG16 . 

SSQ ID MO: 3 gi | 1742273 j dbj | D90778 . 1 | E.coli genomic DNA, Kohara clone #. 

SSQ ID NO: 4 gi | 26107941 | gb j AE016760 . 1 | Escherichia coli CFT073 section . 

SSQ ID NO: 5 gi|41829lemb|X62680.l|ECIS2IS30 E.coli insertion sequences . 

SEQ ID NO: 6 gil24431548|gblAC125613.2| Homo sapiens 3 BAC RP11-379M20 (. 

SSQ ID NO: 7 gij 18159853 |gb|AE009803 . 1 [AE009803 Pyrobaculum aerophilum s. 

SSQ ID NO: 8 gi \ 6067150 [ gb | AC008015 . 5 | AC00801 5 Homo sapiens chromosome 1. 

SSQ ID NO: 9 gi [21644 686 |dbj [AP004363 .3 | Oryza sativa (japonica cultivar. 

SSQ ID NO: 10 gi | 10862756 | emb| AL354852 . 12 | Human DNA sequence from clone 



Value 



Score 


E 


1982 


0.0 


1982 


0.0 


1191 


0.0 


198 


le-47 


44 


0.50 


42 


2.0 


42 


2.0 


42 


2.0 


42 


2.0 



TABLE 2 



{00521 (0053] At this point, those sequences (less than or equal to query-length) that show 
homology/identity with other organisms below a threshold can be identified as unique to the set for which 
they were searched. Such sequences have utility to confirm the presence of at least one member of the 
set, primarily, but not exclusively in a Bioinformatic setting. Sequences unique to higher level sets are 
identified by searching for commonalities between sequences within a classification. These common 
sequences are then searched against the appropriate database. Such sequences, with optional annotation, 
can be realigned with the genomic data, annotated with the information gained saved to a database of 
unique genomic sequences 132, or added to the growing knowledge base of the genome of the organism 
under investigation. 

[00531 C00541 In some embodiments, where a certain minimum number of total unique bases , e.g., 
10,000, is desired across all unique genomic sequences, and less than this number are determined to be 
unique, the genomic data can be preprocessed again, this time dividing the genomic data into smaller 
query-length sequences. The smaller query length sequences are then searched against the target 
database. Control line 133 indicates this path. 

£00^[0055]_The output of the similarity search program can also be used to identify further query- 
length sequences or candidate sequences for organism(s) other than the organism(s) under investigation. 
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For example a first query-length sequence may show high homology/identity only against the particular 
strain it was derived from. But the query sequence might also show some homology to a related strain. 
Such sequences can be referred to as inferred sequences 134. The portion of the related strain where 
limited homology is detected can be searched 120 as a query-length genomic sequence 115 (by being 
searched against the selected database 125) to confirm its identity as a unique genomic sequence 132 for 
the related organism(s); or it can be treated as a candidate sequence for evaluation (discussed below) 
where target-length oligonucleotides are evaluated for amenability to hybridization. Exemplary inferred 
sequences have sufficient homology to the related sequence to be indicated by a BLAST search, but not 
sufficient to cross-hybridize with oligonucleotides derived from the related sequence. 

fQ05»ir00561 Referring to Table 3, a search against the NCBI nr database, using as a query a Vaccinia 
sequence found to be unique by a method of the present invention, identified two candidate sequences 
A01, A02 (a Vaccina strain and the complete Vaccinia genome) with 100% identity over the entire query 
sequence; five non-Vaccinia sequences A03 - A07 with identity ranging from 96% to 92% over portions 
of the query sequence; one non- Vaccinia sequence A08 with 100 identity over a small portions of the 
query sequence; and at least seven non-microbial sequences A09 - A15 with E values greater than three 
for short portions of the query sequence. 

fOO544r00571 T he first group confirms that the query sequence is part of both the Vaccinia strain and 
complete genome. The second and third groups identify sets of organisms with significant homology, 
e.g., E value less than le-5, to the set-unique Vaccinia sequence. Preferred embodiments of the invention 
infer that the second and third group of sequences come from unique regions of the genome of those 
organism sets. Such inferred sequences preferably undergo evaluation and validation as described herein. 
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BLASTN 2.2.4; RID: 1036169670-05727-22152 
Query= (160 letters) 

Database: All GenBank+EMBL+DDBJ+PDB sequences (but no EST, STS, GSS, or phase 0, 1 or 2 
HTGS sequences) 1,430,422 sequences; 7,041,770,514 total letters 
Distribution of 19 BLAST Hits on the Query Sequence 



Sequences producing significant alignments: 



SEP ID NO 



SEP ID NO 



SEP ID NO 



SEP ID NO 



SEP ID NO 



SEP ID NO 



SEP ID NO 



SEP ID NO 



SEP ID NO 



SEP ID NO 



SEP ID NO 



SEP ID NO 



SEP ID NP 



SEP ID NO; 



SEP ID NO 



11 gi!2772662|gb|U94848, 1 IU94848 Vaccinia virus strain Ankara,... 

12 gi|335317|gb|M35027.1|VACCG Vaccinia virus, complete genome 
Hgi|3096962|emb|Yl 1 842. 1 |CVGRI90 Cowpox virus strain GRI-90 
i4_gi|3097015|emb|Y15O35.1|CVY15035 Cowpox virus strain GRI-90... 

15 g i[20152989|gblAF482758.1| Cowpox virus strain Brighton Red... 

16 gij 1 84829 1 3 |gb| AF43 8 1 65 . 1 1 Camelpox virus M-96 from Kazakhs... 

17 gij 1 97 1 7929|gb| AY009089. 1 1 Camelpox virus CMS, complete genome 

18 gi|22123748|gb[AF012825.2[ Ectromelia virus strain Moscow, ... 
i9_gi|14574206|gb|U23449.2| Caenorhabditis elegans cosmid K06A... 
20 g ij687828[gb|U2 13 18.11 Caenorhabditis elegans cosmid K03H9,... 



2J_gi|12000447|gb|AC084754.14| Homo sapiens 12p BAC RP1 1-874G1... 

22 gi|17534934|ref|NM 062895.11 Cuticulin precursor 

23 gi|18250549|emb|AL627429.8| Human DNA sequence from clone R... 
24_gi|16973060|emb|AL590 101.9) Human DNA sequence from clone R... 
25_gi|233372971emb|AL7323 17.13| Mouse DNA sequence from clone ... 

TABLE 3 



Score 


E 


(bits) 


Value 


317 


2e-84 


^17 
j i / 


2e-84 


270 


5e-70 


270 


5e-70 


252 


le-64 


228 


2e-57 


220 


4e-55 


80 


9e-13 


38 


3.2 


38 


3.2 


38 


3.2 


38 


3.2 


38 


3.2 


38 


3.2 


38 


3.2 



{0057^0058] Unique and inferred unique genomic sequences can be identified using the method 
described herein for a number of other biological entities including, but not limited to; Anthrax (Bacillus 
anthracis), Botulism (Clostridium botulinum toxin), Brucellosis (Brucella species), Burkholderia mallei 
(glanders), Burkholderia pseudomallei (melioidosis), Chlamydia psittaci (psittacosis), Cholera (Vibrio 
cholerae), Clostridium perfringens (Epsilon toxin), Coxiella burnetii (Q fever), E. coli 0157:H7 
(Escherichia coli), Emerging infectious diseases such as Nipah virus and hantavirus, Food safety threats 
(e.g., Salmonella species, Escherichia coli 0157:H7, Shigella), Francisella tularensis (tularemia), Ricin 
toxin from Ricinus communis (castor beans), Rickettsia prowazekii (typhus fever), Salmonella Typhi 
(typhoid fever), Salmonellosis (Salmonella species), Smallpox (variola major), Staphylococcal 
enterotoxin B, Variola major (smallpox), Viral encephalitis (alphaviruses e.g., Venezuelan equine 
encephalitis, eastern equine encephalitis, western equine encephalitis), Viral hemorrhagic fevers 
(filoviruses e.g., Ebola, Marburg and arenaviruses e.g., Lassa, Machupo), Yersinia pestis (plague). 
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fQQ5&tf00591 Referring again to Figure 1, unique genomic sequences were realigned 130 with the 
genome from which they were generated in order to determine the start and stop point of each fragment 
within the genome. In preferred embodiments, annotations within the genome in the region containing 
the unique sequence was transferred to the unique sequence. Annotated regions include sequences known 
to have a specific biological function such as protein coding regions, biologically active RNA encoding 
regions, promoter and regulatory elements, spacing elements within operons, protein binding sites, etc. 

f0O#94r00601 In some embodiments of the present invention, the process of obtaining genomic data, 
preprocessing the data, querying the selected database(s) and parsing results to identify candidate 
genomic sequences is implemented as a computer program product. In these embodiments, a plurality of 
organisms and sets of organisms can be investigated concurrently. Computer program products of this 
invention include the ability to indicate the organism(s)/set of organisms of interest, indicate the selected 
database, set thresholds for identifying inferred sequences, direct the handling for inferred sequences, set 
thresholds for identifying unique sequences, direct the handling for unique sequences, output unique 
sequences for evaluation for oligonucleotides amenable to hybridization. Intermediate and final results 
can be made available for user inspection. In preferred embodiments, such computer program products 
are in network communication, e.g., via the Internet with selected databases and databases of genomic 
data such as those available through NCBI via http://www.ncbi.nlm.nih.gov . Operator interface to such 
computer program products is preferably provided through graphical user interface (GUI) technologies as 
known to those skilled in the art. Such computer program products can be configured to operate in a 
centralized fashion or may be distributed over platforms on a network. Some embodiments are 
configured as an application service provider (ASP) accessible by network devices through a browser 

fQQSffpoeil Both unique genomic sequences 132 and select inferred unique sequences 136 are 
evaluated 140 for subsets e.g., target-length oligonucleotides, that are amenable to hybridization. The 
evaluation is done in a target-length oligonucleotide window derived from the query length sequence, and 
preferably moved one base at a time through the query-length sequence. Figure 2 presents a notional 
representation 200 of the incremental progression of a target-length oligonucleotide window 235 through 
a unique genomic sequence 232 one nucletotide 236 at a time. In preferred embodiments, the increment 
of progression is one nucleotide at a time, but in other embodiments different constant or variable 
progressions can be used to take advantage of the inherent properties of the unique genomic sequence. 
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Target-length oligonucleotides are evaluated for, among other characteristics, GC content, metlting 
temperature (T m ), repetitive elements, availability of primer amplification sites, avoiding secondary 
structures such as hairpins and duplexes. In some embodiments this functionality is provided using a 
program such as OLIGO 6 (Molecular Biology Insights, Inc., Cascade CO). In other embodiments, this 
functionality is incorporated into a computer program product of the invention. 

fQQ64¥00621 OLIGO is a multi-functional program that searches for and selects oligonucleotides from 
a sequence file for polymerase chain reaction (PCR), DNA sequencing, site-directed mutagenesis, and 
various hybridization applications. It calculates hybridization temperature and secondary structure of 
oligonucleotides based on the nearest neighbor thermodynamic values. It is also a good tool for 
construction of synthetic genes, finding an appropriate sequencing primer among those already 
synthesized, finding and multiplexing consensus primers and probes, and even finding potential 
restriction sites in a protein. 

ffiO6ar00631 In preferred embodiments, oligonucleotides of approximately 50 bases in length are 
derived from the candidate sequences. A preferred range for oligonucletide lengths is 25 - 100 with a 
range of 50 - 70 being more preferred. Factors that go into determining a range and a preferred value 
include the ability to synthesize the oligonucletoide, the desired hybridization temperature of the 
microarray, balancing melting temperature of the various oligonucleotide against the GC content of the 
molecule and the possible chemical composition of the hybridization solution used on the microarray. As 
these factors change, the preferred length of oligonucleotides will also change. In some embodiments, 
target-length oligonucleotides are chosen based on their melting temperature T m of 90 C, 3'-dimer AG of 
-8.0 kcal/mol, 3 '-terminal stability range of ^1.8 to 11.6 kcal/mol, GC clamp stability of -8.0 kcal/mol, 
minimal acceptable loop AG of -1.9 kcal/mol, maximum number of acceptable sequence repeats of 6 and 
a maximum length of acceptable dimmers of 2 base pairs. In practice, when prepared, oligonucleotides 
are dried and re-hydrated in 3X SSC (a solution of sodium chloride and sodium citrate) at a concentration 
of 1 50 ng/ul. These particular values are the defaults for OLIGO 6 calculations. They can be adjusted by 
the user based on the physical biochemistry of the particular acids. These are good general values for 50- 
mers at this temp and CG content. 
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[00631 TO0641 In some embodiments, favorably evaluated target-length oligonucleotides 145, e.g., those 
found amenable to hybridization, are used as a query to a similarity search program 150, e.g., BLAST. 
The query is directed to a selected database 155 of genome data in order to determine whether the target- 
length oligonucleotide is unique to the organism or organism set under investigation. To this end, 
preferred embodiments parse 150 the similarity search program output to identify oligonucletoides 
lacking significant similarity with other organisms in the selected database, e.g., unique target-length 
oligonucleotides 152. This is counter to the typical use of such search programs. For the purpose of this 
disclosure, "unique" or "uniqueness" is a function of thresholds, preferably controlled by the user, 
regarding identity, homology, score, expected (E) value and the length of the unique sequence under 
consideration. Identity, score, expected value are data returned in a typical BLAST search. In some 
embodiments, lacking significant similarity, e.g., "unique," means no BLAST hits or hits with a E-value 
less than le-5. In some embodiments, this determination of oligonucleotide uniqueness is conducted 
prior to evaluating the oligonucleotides for amenability to hybridization. 

fO0€44r00651 At this point, the certain oligonucleotides that can be identified as unique to the set for 
which they were searched. Such oligonucleotides have utility in confirming the presence of at least one 
member of the set in both Bioinformatic and "wet" settings. Table 4 lists exemplary set-unique 
oligonucleotides for E. coli K12 identified by a method of this invention. Table 5 is an exemplary 
BLAST search result showing the sufficient uniqueness of an oligonucleotide of £. coli 0157:h7. These 
oligonucletides, with optional annotation, can be saved to a database C38 of unique genomic sequence, or 
otherwise added to the growing knowledge base of the genome of the organism under investigation. 

fO0654r0Q661 Referring again to Figure 1, unique oligonucleotide sequences were realigned 160 with 
the genome from which they were generated in order to determine the start and stop point of each 
sequence within the genome. In preferred embodiments, annotations within the genome in the region 
containing the unique sequence was transferred to the unique sequence. Annotated regions include 
sequences known to have a specific biological function such as protein coding regions, biologically active 
RNA encoding regions, promoter and regulatory elements, spacing elements within operons, protein 
binding sites, etc. 
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SEQ ID 

NO: 



26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 



Oligo 
Number 


Organism 


Sequence (5 to 3 ) 




E. coli K12 






NC 000912 


> 29 part354 


91 




ArAnnATATAnTTATArrA^rr,TTATTnTCnTTAGTGGTGAGTCGTCGGT 


92 




rAGTGATAATAACGTGArnrTGGATGGAAAGTTAACTGTTGTATCAGACA 


93 




TATAATrAGAArGGTGAGrTTGTTATTArGGGTGATAATGCCACAGTAAA 


94 




AGAGGTAGAAArrAGCArAnGTGAAAAnAAAAATArrAATGTGACTTGTC_ 


95 




TATrArAGAGnTAGAAArrArrACAGGTnAAAAGAAAAATACCAATGTGA 


96 




TTATArGAnrnTTATTGTrGTTAnTnnTGAnTGnTrGGTATATCTGAATG 


97 




TGAATATGACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAAT 


98 




rAATAGGGTTAATATGAATGGTGGACTTGAAGTGATTGGAGAGAAAAACG 


99 




TTArGGGAGTnATAATAArGTGAGGGTGGATGGAAAGTTAACTGTTGTAT 


100 




ATATTrGrrAnAAGGATTTATTATnrAAAGnAAAACTCTATTGTCGGGCT 


101 




GrAGATATrArAGAGGTAGAAAGCArrArAGGTGAAAAGAAAAATACCAA 


1 02 




AGrGTTATTGTrnTTAGTGGTGAGTrGTrGGTATATrTGAATGGAGATAC 


103 




AAGAGPTTAAArGrTTArrrTrTGAATGrTrTGrTTTAGTCGAACAAAAC 


104 




nAAGTGGCAATATAATrArAArGGTGAGGTTGTTATTACGGGTGATAATG 


105 




rrGGGTTTrTTTATATTrnrrAGAAGGATTTATTATGCAAAGGAAAACTC 



TABLE 4 



BLASTN 2.2.7 [Jan-02-2004] 

RID: 1074620345-32204-105313520645.BLASTQ4 
Query= (50 letters) 

Database: All GenBank+EMBL+DDBJ+PDB sequences (but no EST, STS, GSS, or phase 0, 1 or 2 HTGS sequences) 2,017,250 sequences; 
9,771,1 19,756 total letters 

Distribution of 1 10 Blast Hits on the Query Sequence 

Score E 

Sequences producing significant alignments: (bits) 

SEP ID NO: 41 g ill2519298igb|AE005660.IlAE005660 Escherichia coli 0157:H7... 100 

SEP ID NO: 42 gill3364704ldbilAP002569.il Escherichia coli 0157:H7 DNA, c... 100 
SEP ID NO: 43 g i|262091|gb|S50878.1|S50878 tyrosine hydroxylase {5Yegion... 

SEP ID NO: 44 g i[244302661emblAL928973.3| Mouse DNA sequence from clone R... 38 

SEP ID NO: 45 g i|38083308{ref]XM 1 1 1 684.3| Mus musculus similar to corneo... 36 

SEP JD NO: 46 g i|380831 12|ref1XM 359287.11 Mus musculus similar to comeo... 36 

SEP ID NP: 47 gi l332869071gblBC055373.1| Mus musculus cDNA clone MGO6076... 36 
SEP ID NO: 48 g i|18412479|reflNM 106620.11 Arabidopsis thaliana hypothetic. 



TABLE 5 

£006gI00621_Just as set-unique genomic sequences and oligonucleotides can be identified for strains of 
a particular organisms (E. coli as above), the method can be used to identify genomic sequences and 
oligonucleotides unique to higher level sets such as species, genus, family, clade, ad hoc sets, etc. 



40 



36 



3e-19 
3e-19 

0.95 
3.8 
3.8 
3.8 



0.24 



3.8 
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fOO ^r00681 Embodiments of the invention include both a method for profiling the hybridization 
response of set-unique target oligonucleotide sequences to a variety of organisms and sets of organisms 
and the resulting database of profiles. In these embodiments, target oligonucleotides are created in 
accordance with the unique sequences and are spotted to a microarray. Genomic DNA from each 
organism of interest is purified using procedures known to those skilled in the art. Preferably, the 
genomic DNA includes DNA from the area surrounding the region from which at least one target-length 
oligonucleotide was derived. In preferred embodiments, aliquots of the DNA are labeled, preferably with 
fluorescent dNTPs in a Klenow reaction. Each labeled DNA sample is allowed to hybridize with a 
microarray of the target-length oligonucleotides. The hybridized microarrays are wash, scanned, and the 
data is imported into a data visualization program (e.g., the suite of analysis software offered by Spotfire 
of Somerville MA). The data is evaluated to determine that target-length oligonucleotides exhibit 
true_positive and truenegative reactions to each organism of interest. 

fQ96SP0691 In practice, multiple copies of target-length oligonucleotides and multiple trials of 
hybridization are used for each organism of interest. In this fashion, data can be collected on the variation 
of hybridization intensity to be expected between each target-length oligonucleotide and each organism of 
interest. For example, statistics regarding the distribution of hybridization intensity between a target- 
length oligonucleotide and its corresponding organism can be collected. Table 6 presents a notional 
example of hybridization intensity for two organisms against oligonucleotides found to be unique to 
separate strains of E. coli and the E. coli Shiga gene. Table 6 identifies exemplary unique oligonucleotide 
sequences and provides the hybridization intensity for those sequences showing the greatest differential in 
hybridization intensity as between E. coli K 12 and E. coli ol57:h7 based on the query-length sequences 
NC_009 13_29 _part354 and NC_002695_194_part29 respectively, along with Shiga gene oligonucleotide 
that exhibit differential hybridization intensity between the two £. coli strains. Figure 3 illustrates this 
differential hybridization intensity between strains for oligonucleotides indicated by parenthetical 
numbers in Table 6. 
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SEQ 



Oligo No. 


Organism 


Sequence (5' to 3') 


K12 

Intensity 


0157:H7 
Intensity 




E. coli K12, NC 


000913 29 part354 






(l) 91 




ACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGT 


33298 


13380 


92 




CAGT G AT AAT AACG TG ACG CT GG ATGG AAAGTT AACTG TTGT ATCAG ACA 






(2) 93 




TATAATCACAACXX5TGAGCTTGTTATTACGGGTGATAATGCCACAGTAAA 


13973 


5880 


94 




AGAG^TAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTC 






95 




TATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGA 






(3) 96 




TTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATG 


57432 


25699 


97 




TGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAAT 






(4) 98 




CAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACG 


16711 


6814 


99 




TTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTAT 






(5) 100 




ATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCT 


17384 


5462 


102 




AGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATAC 






(6) 104 




GAAGTGGCAATATAATCAC^CGCTGAGCTTGTTATTACGGGTGATAATG 


13862 


2643J 


105 




CCGGCI^TC^TATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTC 








E.coli 0157:h7 


, NC 002695 194 part29 






278 








| 


279 




ATTGTAACX5CGATTACTTTCTCAATCTTCTGAAATACGCAGTGCAAGTAG 






281 




TGTAAGGCGATTACTTTCTCAATCTTCTGAAATACGCAGTGCAAGTAGTC 






283 




TTAGGTCACACAATTCCJWVAGGCX^TTATGAACTTACTGAAGAACTTCC 






285 




CJ^TCTTCTGAAATACGCAGTGCAAGTAGTCTTCAGGTTGAAAGAATCTG 






(7) 287 




GTTATGAACTTACTGAAGAACTTCCCTGGGTTTTTGAAGCGATGAATGAG 


6038 


14850 


288 




GTCO^TTCTGTTAGGTCACACAATTCCAAAAGGCGGTTATGAACTTACT 






(8) 289 




GTAGTCTTCAGGTTGAAAGAATCTGGGTCGATGTTGATGACAGTTGGTAT 


14399 


4/ I 


292 




AAATCTCATTGTCX^CTTCATGAAGGTCCAATTCTGT^AGGTCACACAATT 








E. coli Shiga G 


ene AF461172 






296 










297 




CTGACrATCATGGACAAGACTCTGTTCGTGTAGGAAGAAi i iVn -iTGGA 






(9) 298 




AGGTACAACAGCX^TTACATTGTCTGGTGACAGTAGCTATACCACGTTAC 


9529 


18509 


299 




GTGAGCTATACGGAAAGTACACAAAAAGGAAGGTGCGACCACAATTAATA 


17781 


27191 


) 300 




ATr^CC^TTCGTTGACTACTTCTTATCTGGATTTAATGTCGCATAGTGGA 






313 




CGG AAAGT A CA CAAAAAGGAAGGTG CGA CCACAAT T AATAACAAAAT CTT 


7644 


20573 


1 314 




CGCTCTGCJ^TAGGTACTCCATTACAGACTATTTCATCAGGAGGTACGTC 






\ 322 




ATTTACCAACAGATGGAATCTTCAGTCTCTTCTTCTCAGTGCGCAAATTA 






1 323 




AAGGAAGGTGCGACCACAATTAATAACAAAATCTTAAAAATTGCACATGG 






) 324 




AATTATTTACCAACAGATGGAATCTTCAGTCTCTTCTTCTCAGTGCGCAA 


5319 


20729 


» 325 




AAGTTATTTTTCGTTGACTCAGAATAGCTCAGTGAAAATAGCAGGCGGAG 


13167 


37404 



TABLE 6 
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[00 6 91 f 00 701 In one embodiment, the genomes of both E. coli K12 and E. coli 0157:H7 were 

investigated to identify unique sequence. This search included queries such as the query-length sequence 
in Tabid. Sequences such as NC_000913_29_part354 (Table l)and NC_002695_194j)art29 and the E. 
coli Shiga Gene AF461 172 were BLAST searched against the entire NCBI nt database to verify that these 
sequences are unique as defined with respect to embodiments of the invention, e.g„ see Table 2 for the E. 
coli K12 results. These results confirmed that the three query-length sequences noted directly above are 
unique. The query length sequences were used to generate oligonucleotides. A exemplary list of these 
oligonucleotides is presented in Table 4 for E. coli K12. The full set of oligonucleotides were BLAST 
searched against the NCBI nr database to verify that they were unique at the oligonucleotide level; Table 
5 presents an exemplary extract of the BLAST output for an oligonucleotide of E. coli ol 57:57. 
Oligonucleotides confirmed unique for the set were manufactured and spotted to microarrays. The 
microarrays were hybridized with genomic probes manufactured by Klenow labeling genomic DNA with 
cy3-dCTP. Oligonucleotides that demonstrated differential hybridization patterns were detected; see 
Table 6 for exemplary values and Figure 3 for a graphic representation. The thirteen oligonucleotides can 
now be used to distinguish between E. coli K12 and E. coli 0157:57. 

r0070i r00711 In some embodiments, the invention leverages the principles illustrated by the previous 
example. In those embodiments, the database of hybridization profiles is used as one source to pick 
oligonucleotides that are informative at decision points in decision tree. Structures appropriate for the 
decision tree include most taxonomic hierarchies, but any hierarchy where some oligonucleotides at 
sibling decision points have discernable differential hybridization will work. Placement of redundant 
unique and conserved phylogenetic specific oligonucleotide targets permits the identification of sets of 
organisms (e.g., family, class, order, genus, species, strain) specific branch points that can be used to 
identify organisms. In addition, placement of targets can include those sensitive to other features such as 
virulence genes, structural genes, biochemical genes, antibiotic resistance genes, housekeeping genes, 

f0O74 : }f0072] Figure 4 illustrates a phylogenetic decision tree 400 in accordance with embodiments of 
the invention. In this particular embodiment, multiple oligonucleotide targets (including those from a 
database of unique oligonucleotides as described above that were found to provide useful differential 
hybridization between the organism sets under consideration) were spotted onto a microarray for each 
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decision pint, e.g., 401 and exposed to a labeled complex sample for hybridization. Each decision point, 
e.g., 401 is scored with absent (-) or present (+) where: 

Score = ±* 

T 

A = number of (+) hybridized data points; 
B = number of (-) hybridized data points; and 
T = total data points. 

feQ7«r00731 In the embodiment illustrated in Figure 4, decision points having a score above a 
threshold are indicated as present (+). In other embodiments, actual hybridization values are reported. In 
further embodiments, absent/present calls are either made on the basis of statistics for all hybridization 
points associated with the decision point. In yet other embodiments, the node with highest score is called 
"present" under a branch point. 

[00731 IW741 In some embodiments, a confidence score is determined. In general each successive 
present call through the phylogenetic tree increases the confidence of the call below it and the confidence 
of the final identification. In Figure 4, the presence of two present calls for two different strains of E. coli 
decreases the confidence of an accurate identification at the strain level but increases the confidence of 
the identification at the level of the species. In some embodiments, the confidence score for a decision 
point is determined based on Bayesian analysis methodologies where present calls in the correct lineage 
contribute to the confidence score for the final identification. Factors that increase the confidence of the 
final identification include present calls in the phylogenetic positions in the tier above the absent/present 
call and additional absent calls laterally within that tier. Also contributing to the confidence of the final 
identification is the presence of lateral absent calls within the same tier as the absent/present call. Thus 
the absent/present call at any one location in the phylogenetic tree is dependent on all of the 
absent/present calls in all tiers. 

fOO744r00751 I n some embodiments, the confidence score is represented as a numerical percentile 
between 0% and 100% by the following equation 
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Confidence _ score = 



x\00, where 



T 



a = a + c + e, and 



p = b + d + f , where 



a = number of correct (+) calls in the correct lineage 
b = number of incorrect (-) calls in the correct lineage; 
c = number of correct (-) calls in the correct lineage; 
d = number of incorrect (-) calls in the correct lineage; 
e = number of correct (-) calls in the incorrect lineage; and 
/= number of incorrect (+) calls in the incorrect lineage. 

[00751 [00761 In some embodiments, a virulence score, a value used in conjunction with the final 
identification of the organism, is associated with the relative power and degree of pathogenicity of an 
organisms to produce disease. This value is derived from the number of virulence factors present, as 
detected by hybridization, compared to the total number of virulence factors represented in the array for a 
particular organism. Virulence targets can be spotted to the microarray independent of targets for 
absence/presence. This score is represented as a numerical value between -1 and +1 by the following 



[00761 [00771 I n addition, the ability to identify virulence versus structural components by 
oligonucleotide-specific hybridization permits the identification of recombinant organisms that contain 
structural components of one organism and the virulence components of a different organism. In some 
situations, a large number of data points may be required to identify all species and strains of micro 
organisms that might be found in a complex biological sample. 



equation: 



Virulence score - - , where 




r 



the number of (+) hybridized virulence-specific data points, and 

the number of (-) hybridized virulence-specific data points. 

the total number of virulence-specific targets associated with the decision point. 
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fQO7?4r007S1 Some methods of the present invention are conducted to mitigate the effect of "noise" by 
pre-screening. In accordance with these methods, a background sample of interest is obtained, and 
nucleic acid sequences in the sample are amplified using random amplification and combined with a 
microarray for hybridization. Signals from amplification products that hybridize with the microarray are 
recorded to be discounted in subsequent analysis. 

fflQ3CT00791 These methods are, for example, suited to customs applications. Customs officials at 
ports of entry including airports, harbors, and country borders can utilize prescreening to screen food 
samples for commonly occurring pathogens such as E. coli, Salmonella typhi, Hepatitis A virus and the 
like. In pathogen-free samples the level of hybridization observed to know pathogens on the array is 
minimal, the hybridization profile of such pathogen-free samples is then used as a baseline level to 
subsequently identify contaminated samples. 
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ABSTRACT 

Identifying genomic sequences and oligonucleotide sequences unique to a set of organisms. 
Methods include obtaining genomic data characteristic of the set; formatting the data into at least one 
query-length sequence, each query-length sequence being of a format compatible with a similarity search 
enginer-Sisearching a selected genomic database using the query and the similarity search engine^-F: and 
then parsing the results of the search for those sequences showing uniqueness to the set. 
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